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Abstract 

The Subtractive Optimally Localized Averages (SOLA) 
method, developed and extensively used in helioscismol- 
ogy, is applied to artificial data to obtain measures of the 
sound speed inside a solar-type star. In contrast to inver- 
sion methods which fit models to some aspect of the data, 
methods such as SOLA provide an honest assessment of 
what can truly be resolved using seismic data, without in- 
troducing additional assumptions such as that the space 
of admissible stellar models is small. The resulting mea- 
sures obtained from SOLA inversion can subsequently be 
used to eliminate putative stellar models. Here we present 
results of experiments to test the reliability of SOLA in- 
ferences using solar models and models of solar-type stars. 
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1. Introduction 

Seismology of solar-type stars is expected in the not-too- 
distant future to provide information of relevance for un- 
derstanding stellar structure. We need to develop tools 
that will allow us to use the full potential of the informa- 
tion provided by the seismic observations of other stars. 

Early attempts at inverting artificial sets of low-degree 
modes have given mixed results. Gough & Kosovichev 
(1993) and Roxburgh et al. (1998) inverted the frequen- 
cies calculated for a 1.1 M© and a 0.8 M star, respec- 
tively to obtained very encouraging results. More recently, 
Berthomieu et al. (2001) carried out a careful analysis of 
the results to be expected in inversions of solar-like os- 
cillations, taking into account the expected errors in the 
COROT observations and mode amplitudes in a 1.45 M© 
star; they concluded that reasonably well-localised inver- 
sion is possible in the core of such a star. On the other 
hand, Basu et al. (2001) failed in their attempts at invert- 
ing a set of low-degree modes of a solar-type star. 

One of the main differences between these efforts was 
the mode set used for the inversions. The Roxburgh et 
al. work used a very optimistic set: it is unlikely that all 
those modes can be observed. The Gough & Kosovichev 
mode set was somewhat more conservative. Basu et al. 
(2001) used two conservative mode sets. The modes sets 
and errors used in the different works are summarized in 



Table. 1. In addition, unlike Gough & Kosovichev and 
Roxburgh et al., Basu et al. tried to invert for the sound 
speed c using density p as the second variable, instead of 
inverting for u = p/p, p being the pressure, with the he- 
lium abundance Y as the second variable. It is known that 
because density kernels (at fixed c 2 ) have much larger am- 
plitudes than Y kernels (at fixed u) , a c inversion is more 
difficult than a u inversion. The Basu et al. (2001) inver- 
sions were done under the assumption that the mass, or 
radius, or both of these properties of the test star was 
not known very well and hence this introduced additional 
uncertainty. Also we believe that there were numerical er- 
rors in the kernels used. 

In this work we attempt to understand the differences 
in the results, in particular we study the effect of the mode 
set and errors on the results obtained. The mode sets used 
are the same as those used by Basu et al. (2001); the mode 
set of Gough & Kosovichev (henceforth the GK set) is used 
for comparison. In order to be able to compare with older 
results we try to invert for u rather than c 2 which is the 
variable of choice in solar inversions. To avoid additional 
uncertainties due to uncertainty in the radius or the mass 
of the star under consideration, in this work we confine 
our attention to solar models. Finally, we note that the 
assumed errors are much larger than is typical in solar 
frequencies. 

2. Models and inversion techniques 

Inversions for stellar structure are based on linearizing the 
equations of stellar oscillations around a known reference 
model. The relative differences in u (i.e. Su/u) and the 
difference in Y (5Y) between a star and a reference stellar 
model can be related to the differences in the frequencies 
of the star and the model [6uJi/u>i), 

5 -^ = J < y (r)dr + / K{.Jr)Ar + (1) 

where r is the normalized distance to the centre. The in- 
dex i numbers the multiplets (n, I). The kernels and 
K l n 2 are known functions of the reference model. The 
term in F sur f(uji) is the contribution from the uncertain- 
ties in the near-surface region (e.g. Christensen-Dalsgaard 
& Berthomieu 1991); here Qi is the mode inertia, normal- 
ized by the inertia of a radial mode of the same frequency. 
In general, the right-hand side of equation (1) may also 
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Table 1. Mode sets considered and assumed standard de- 
viation of errors. 



Paper 


Mode set 


Assumed 






errors 


Gough & Kosovichev 


1 = 0,1, n=10-30 


0.3 uRz 


(GK set) 


/ = 2, n=9-29 




Basu et al. 


1 = 0, n =11-27 


0.1a*Hz 


Set 1 


1 = 1, n =12-28 






1 = 2, n =15-27 




Basu et al. 


1 = 0, n =14-32 


0.3 /jHz* 


Set 2 


1 = 1, n =13-29 






1 = 2, n =15-30 






1 = 3, n =16-28 




Roxburgh et al. 


I = 0, n=l-31 


0.3 /iHz 




1 = 1,2, n =1-30 






1 = 3, n =1-29 





Here we use the Subtractive Optimally Localized Av- 
erages (SOLA) method (Pijpers & Thompson 1992, 1994) 
to determine the inversion coefficients such that the aver- 
aging kernel is an approximation to a given target T(ro, r). 
Details on the application of the SOLA technique to struc- 
ture inversion were given by, e.g., Rabello-Soares, Basu & 
Christensen-Dalsgaard (1999). 

The reference model used in this work is the solar 
Model S of Christensen-Dalsgaard et al. (1996). This is 
a standard solar model. 

The test model (the proxy star) we use is model MIX 
of Basu, Pinsonneault & Bahcall (2000). This is a non- 
standard model, with an artificially mixed core. This model 
was selected because the differences in u with respect to 
the reference model are large, somewhat along the lines 
of what we expect for stars other than the Sun due to 
uncertainties in their ages. 



* For Set 2, various errors are considered (see Fig. 5). 

contain a term ^S\n(M/R 3 ) to absorb scaling of the fre- 
quencies with stellar mass M and radius R according to 
(GM/R 3 ) 1 / 2 . 

For linear inversion methods, the solution at a given 
point ro is determined by a set of inversion coefficients 
Ci( r o), such that the inferred value of, say, 5u/u is 

* ' i 

From the corresponding linear combination of equations 
it follows that the solution is characterized by the av- 
eraging kernel, obtained as 

K(r ,r) = J2*(ro)Ki Y (r), (3) 

i 

and also by the cross-term kernel: 

C(r ,r)=Y, c ^ K YJr) , (4) 

i 

which measures the influence of the contribution from SY 
on the inferred du/u. 

The surface term in equation ([!]) may be suppressed 
by assuming that i^urf can be expanded in terms of poly- 
nomials tpxt an d constraining the inversion coefficients to 
satisfy £\ Ci(r Q )QT x il)\(Ui) = 0, A = 0,1,..., A. (Dappen 
et al. 1991). 

The goal of the inversions is to obtain a localized aver- 
aging kernel, while suppressing the contributions from the 
cross term and the surface term in the linear combination 
in equation (0), and limiting the error in the solution. 
Also, /C(ro, r) must have unit integral with respect to r. If 
this can be achieved, then 

(^(r ))./x;(ro,r)^dr (5) 
defines a proper average of du/u. 



3. Results 



30 



. Comparing c 2 and u inversions 


. ' ' ' i ' ' ' i ' ' ' i ' ' ' i ' ' ' 

"- Averaging kernel for c 2 

: A Averaging kernel for u 




? \ (a) r„ = 0.05 R. 








: 1 1 1 -f- 1 1 1 1 1 1 1 1 1 1 1 1 
1- (b) r =0.139 R, 




EH— 1 — H-^H — 1 — 1— | — 1— 1 — 1 — — 1 — 1 — 1— | — 1— 1 — h- 




L (c) r = 0.227 R. 









0.2 0.4 0.6 0.8 1 



Figure 1. Averaging kernels for sound-speed inversions (con- 
tinuous) and u inversions (dashed) for a few target radii. The 
mode set used was Set 2 with a uniform error of 0.3 uHz. 

Fig. ^ shows a comparison of averaging kernels ob- 
tained for c 2 and u inversions. The data set used was Set 2, 
with uniform errors of 0.3 /iHz. Note that the c 2 averaging 
kernels have a lot of surface structure, which has a detri- 
mental effect on the inversion results. The mode set needs 
to be expanded in order to get a cleaner averaging kernel. 
Fig. H shows a comparison of the cross-term kernels ob- 
tained for c 2 and u inversions. One can see that with the 
mode set used (Set 2), we cannot successfully suppress the 
cross term in a c 2 inversion. In contrast, the cross term for 
the u inversions is very small. 
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Figure 2. Cross-term kernels for the inversions. r/R, r/R. r/R. 



3.2. Effect of mode set 
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Figure 3. The averaging kernels obtained with Set 1 (dotted 
lines), Set 2 (continuous lines) and GK (dashed lines) mode 
sets at 0.05 (Panel a), 0.139 R* (Panel b) and 0.227 
(Panel c). The near surface structure of these kernels is shown 
in the panels to the right. 

Fig. H shows a comparison of the averaging kernels 
obtained for u inversions using mode sets Set 1, Set 2 and 
GK. A uniform error of 0.3 /iHz was assumed for all sets. 
Note that the Set 1 averaging kernels have structure at 
the surface. Set 2 and GK kernels are much cleaner. Set 1 
also results in very large cross-term kernels. It should be 
noted that MOLA inversions can give somewhat better 
averaging kernels. 

Fig. H shows the inversion results for the three sets. The 
horizontal error bars are an indication of the resolution of 



Figure 4- Results of u inversion for Set 1 (panel a), Set 2 
(panel b) and GK (panel c). The continuous line is the ex- 
act difference between the models. Inversion results are plot- 
ted at the corresponding target radius. The vertical error bars 
show the 1-a uncertainty in the inversion result, obtained by 
propagating the mode uncertainties through the inversion pro- 
cess. The horizontal error bars provide an indication of the 
resolution of the inversion; they extend from the first to the 
third quartile point of the averaging kernels. A uniform error 
of 0.3 fiHz was assumed for each set. 

the inversion. One can see that the Set 1 results have very 
large errors, as well as poor resolution. Set 2 and GK give 
very good results for a substantial portion of the core. 

3.3. Effect of mode errors 

A reduction in mode errors helps the inversions. The ef- 
fects can be felt in two ways: one could improve the res- 
olution, while keeping the error on the solution the same, 
or, one could keep the averaging kernel the same and de- 
crease the error on the solution. Fig. || shows the inversion 
results for Set 2 for three different error-distributions: (1) 
uniform errors of 0.3 /iHz, (2) uniform errors of 0.1 /iHz, 
and (3) an error-distribution which mimics the distribu- 
tion of errors in solar oscillations frequencies. The distri- 
bution was scaled to be around 0.1 /iHz at low frequency 
(see Fig. || for the error distribution). The reduction of 
errors gives a more dramatic improvement for Set 1 re- 
sults. If the errors are reduced by a factor of 5, the surface 
structure of the averaging kernels can be removed and the 
cross-term kernels can be made much smaller. 

4. Discussion 

In common with other investigations of the resolving power 
of the low-degree modes expected from observations of 
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Figure 5. Results ofu inversion with Set 2. Different error dis- 
tributions were used for the inversions whose results are shown 
in the different panels. 
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Figure 6. The distribution of errors used for the inversions 
shown in Fig. 5(c). The errors were scaled from the error distri- 
bution of solar oscillation frequencies so that the low-frequency 
part of the distribution is around 0.1 fiHz. 

solar-like oscillations in distant stars, we find that with 
realistic assumptions about the expected mode set and 
observational errors it is possible to achieve some limited 
resolution in the core of the model. With our present proxy 
star this, for example, shows clear evidence for the core 
mixing as reflected in the steep rise in Su/u towards the 
centre. 

The results depend strongly on the choice of variables 
in the inversion. For helioseismic inversion, a common 
choice is to express the inverse problem in terms of adi- 
abatic sound speed and density. With the limited set of 



modes available in the stellar case, the quality of the inver- 
sion for the sound speed suffers from the need to suppress 
the contribution from the density. A preferable choice of 
variables is the pair (u, Y); here the contribution from Y 
is small and essentially confined to the outer layers, where 
helium is partly ionized. 

Not surprisingly, the results show substantial depen- 
dence on the choice of mode set, although the precise 
manner in which this happens is not obvious and needs 
further study. It is important here also to include more 
realistic estimates of the expected properties of the modes 
and observations, as has been attempted by Berthomieu 
ct al. (2001). 

In analyses of evolved stars it is possible that modes 
with partial g-mode character may be observed; some evi- 
dence for this may already have been found in the subgiant 
f] Bootis (e.g. Christensen-Dalsgaard, Bedding & Kjeldscn 
1995). Such modes would have a highly beneficial impact 
on inversions for the stellar core, because the averaging 
kernels could exploit the localized trapping of the eigen- 
functions. 
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